Dynamics of a FitzHugh-Nagumo system subjected to 
autocorrelated noise 
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We analyze the dynamics of the FitzHugh-Nagumo (FHN) model in the presence 
of colored noise and a periodic signal. Two cases are considered: (i) the dynamics 
of the membrane potential is affected by the noise, (ii) the slow dynamics of the 
recovery variable is subject to noise. We investigate the role of the colored noise 
on the neuron dynamics by the mean response time (MKT) of the neuron. We 
find meaningful modifications of the resonant activation (RA) and noise enhanced 
stability (NES) phenomena due to the correlation time of the noise. For strongly 
correlated noise we observe suppression of NES effect and persistence of RA phe- 
nomenon, with an efficiency enhancement of the neuronal response. Finally we 
show that the self-correlation of the colored noise causes a reduction of the effective 
noise intensity, which appears as a rescaling of the fluctuations affecting the FHN 
system. 
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I. INTRODUCTION 



Physical and biological systems are con- 
tinuously perturbed by random fluctuations 
produced by noise sources always present in 
open systems (]l|). In physical and biolog- 
ical systems, noise can be responsible for 
several interesting and counterintuitive ef- 
fects. Resonant activation (RA) and noise 
enhanced stability (NES) are two noise ac- 
tivated phenomena that have been widely 
investigated in a dissipative optical lattice, 
spin systems, Josephson junctions, chemical 
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and biological complex systems, and financial 
markets |1; H; 0; [s^. Because of the stochas- 
tic resonance phenomenon, a weak electric 
stimulus can be enhanced up to the detec- 
tion threshold of the neural system (jl). 

Recently the increasing interest in neu- 
ronal dynamics gave rise to many studies in 
this field. Two approaches, widely used to in- 
vestigate the neuron response to an external 
stimulus, are the Ho dgkin- Huxley model (0) 
and its simplified version, the FitzHugh- 
Nagumo (FHN) model (0; i; M, [HI; E B 
Q). A FHN system, in the presence of fluc- 
tuations, because of its intrinsic nonlinearity, 
can give rise to interesting noise induced phe- 
nomena: modification of detection threshold 
by manipulation of noisy parameters (jisl ). 
noise-induced activation and coherence reso- 



nance (jl6l ). stochastic resonance in the pres- 
ence of colored noise with a 1//^ spectrum 
(0 < /3 < 2) (jirl). optimization of the inter- 
spike time and dependence of neuron firing 
on subthreshold periodic signal and additive 
noise (1181). resonant activation and noise en- 



hanced stability (jl9l ). 

Despite the broad interest in FHN system, 
however, understanding the role of the fluctu- 
ations and their correlation in nerve cells still 
remain elusive. This is why it is of particu- 
lar importance the investigation of the mean 
response time of a single neuron-like element 
function of the noise parameters in the 
presence of various external perturbations. 
In this paper, by using a colored noise given 
by the archetypal Ornstein-Uhlenbeck pro- 
cess, with correlation time r, we study the 
role of the fluctuations on the FHN system. 
In particular we analyze the mean response 
time (MRT) as a function both of the noise 
intensity and the correlation time and we 
find that the RA and NES phenomena ob- 
served in the presence of white noise ( llOl ) 
show meaningful modifications when a col- 
ored noise source is used. The paper is or- 
ganized as follows. In the first paragraph we 
shortly introduce the Ho dgkin- Huxley (HH) 
system, the archetypal model to describe neu- 
ronal dynamics, and the FitzHugh-Nagumo 
(FHN) system, which is a simplified version 
of the HH model, discussing briefly the mo- 
tivations for using, instead of the more real- 
istic HH approach, the FHN model. In order 
to analyze the dynamics of a single neuron, 
we take into account a FHN system driven 
by a periodic signal and we analyze the de- 
terministic dynamics. Afterwards we investi- 
gate the stochastic FHN system. In particu- 
lar we consider the effects of the fluctuations 
by inserting in the deterministic FHN model 
a colored noise source. We analyze two dif- 
ferent cases: the noise affects the dynamics 
of the membrane potential (case I), the re- 
covery variable is subject to noise (case II). 
We introduce the mean response time (MRT) 
in order to calculate the response time of the 



neuron, that is the time that the membrane 
potential takes to reach a given threshold. Fi- 
nally we show that the self-correlation of the 
colored noise causes a reduction of the effec- 
tive noise intensity, that is a rescaling effect 
of the fluctuations affecting the FHN system. 
In the last section we draw the conclusions. 



II. HODGKIN-HUXLEY MODEL 

In the year 1952 Alan Lloyd Hodgkin and 
Andrew Huxley proposed a neuronal model 
to describe the dynamics of an action poten- 
tial. This consists in a "spike" of electrical 
discharge generated as response to an exter- 
nal stimulus and travelling along the mem- 
brane of a cell. The Ho dgkin- Huxley (HH) 
model, which consists of four nonlinear ordi- 
nary differential equations, describes the elec- 
trical characteristics of excitable cells such as 
neurons and cardiac myocytes, explaining the 
ionic mechanisms responsible for the appear- 
ance and propagation of action potentials in 
the squid giant axon (7). The semiperme- 
able cell membrane separates the interior of 
the cell from the extracellular liquid and acts 
ClS db capacitor. If an input current I{t) is in- 
jected into the cell, it may add further charge 
on the capacitor, or leaks through the chan- 
nels in the cell membrane. Because of ac- 
tive ion transport through the cell membrane, 
the ion concentration inside the cell is dif- 
ferent from that in the extracellular liquid. 
The Nernst potential generated by the differ- 
ence in ion concentration is represented by a 
battery. Therefore the Hodgkin-Huxley type 
model represents the biophysical characteris- 
tic of cell membranes. 

The HH model consists of four equations. 
The first one, given by 
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du 



(1) 



with 



^h = gNam^h{u-ENa)+gKn^{u-EK)+gL{u-EL] 

k 

(2) 
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provides the dynamics of the vohage u across 
the membrane and Ik is the sum of the ionic 
currents which pass through the cell mem- 
brane. The conservation of electric charge on 
a piece of membrane implies that the applied 
current J(t) may be split in a capacitive cur- 
rent Ic = Cdu/dt which charges the capaci- 
tor C and further components Ik which pass 
through the ion channels: a sodium channel 
with index Na, a potassium channel with in- 
dex K and an unspecific leakage channel with 
resistance R. All channels may be character- 
ized by their resistance or, equivalently, by 
their conductance. The leakage channel is 
described by a voltage-independent conduc- 
tance ql = i/R] the conductance of the other 
ion channels is voltage and time dependent. 
If all channels are open, they transmit cur- 
rents with a maximum conductance QNa or 
gx, respectively. The parameters E^^a, Ex, 
and El are the reversal potentials. Rever- 
sal potentials and conductances are empirical 
parameters whose values, referred to a volt- 
age scale where the resting potential is zero, 
were obtained by Hodgkin and Huxley (0). 
The three additional variables m, n, and h 
describe the probabilities that the two chan- 
nels are open. The combined action of m and 
h controls the Na'^ channels. The gates 
are controlled by n. The three variables m, 
n, and h are called gating variables. They 
evolve according to the differential equations 

d/TTl 

— = am{u){l - m) - I3m{u)m (3) 
at 

— = - n) - f3n{u)n (4) 
dt 

dh 

— = -h)~ ^u{u)h. (5) 

The functions a.^, On, (^h, Pm, Pn, Ph are em- 
pirical functions of u that have been adjusted 
by Hodgkin and Huxley to fit the data of 
the giant axon of the squid. The four equa- 
tions ([I]),© -([5]) define the Hodgkin-Huxley 
model. 



III. FITZHUGH-NAGUMO MODEL 

The FitzHugh-Nagumo (FHN) model is 
a simpler version of the Hodgkin-Huxley 
model. In particular the FHN model takes 
into account the excitable variable, that is 
the membrane potential, which exhibits a fast 
dynamics, and the recovery variable, char- 
acterized by a slow dynamics and responsi- 
ble for the refractory behaviour of the neu- 
ron. The reduction from four variables (HH 
model) to two variables (FHN model) allows 
to focalize the analysis on the properties of 
excitation and propagation, neglecting the 
specific electrochemical properties of sodium 
and potassium ion flow. Thus the neuron dy- 
namics can be described in terms of a voltage- 
like variable and a recovery variable. The for- 
mer, governed by a cubic nonlinearity, pro- 
vides a regenerative self-excitation by a pos- 
itive feedback, the latter, subject to a linear 
dynamics, is responsible for a slower negative 
feedback. 

The archetypal form for the FHN model 
is given by the following two equations 

V = f{V) -W + I (6) 
W = a{bV - cW) + d (7) 

where f{V) is a polynomial of third degree, 
and a, b, c, d are parameters with constant 
values. 

Actually the Hodgkin-Huxley model is closer 
to the real behaviour of neuronal dynam- 
ics. However, only projections of its four- 
dimensional phase trajectories can be ob- 
served. On the other side, by using the 
FitzHugh-Nagumo model, the solution for 
the time behaviour of the neuronal response 
can be represented in a two-dimensional 
space. This permits to give a geometri- 
cal explanation of important biological phe- 
nomena connected with neuronal excitability 
and spike-generating mechanism. Eqs.(l6]),(I7]) 
provide the deterministic dynamics for a 
FHN neuron. In fact, in the presence of an 
external stimulus, with enough strong ampli- 
tude, the membrane voltage (fast variable) 
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FIG. 1 Phase diagram for FitzHugh-Nagumo 
model (modified from FitzHugh 1961 (0)). The 
(V,W) trajectory (solid line) after a spike event 
is in a refractory state, before reaching the equi- 
librium state which corresponds to the point 
where the two nullclines (dash line for V = 
and dash-dot line for W = 0) intersect. Here we 
set b/c = 1.25, d/{ac) = 0.875, / = -0.001. 



increases rapidly, crossing a threshold. This 
event (voltage spiking) represents the neuron 
firing, that is the nerve response due to the 
external stimulus. After spike generation the 
voltage goes rapidly under threshold and a re- 
fractory time occurs so that during this time 
no further firing is possible. This behaviour 
is represented in Fig. [T] where a (V", W) tra- 
jectory and the two nullclines 



W 



W 



fiv)- 



d 



ac 



are shown for f{V) = V — 



(9) 
1.25, 



^ = 0.875, / 



-0.001. 



IV. FHN MODEL WITH PERIODICAL 
DRIVING SIGNAL 

In the presence of a driving signal the FHN 
model becomes 

X = X — - — y + Asinioj t + ifo) (10) 




FIG. 2 Parameter plane {uj,A) for = 0. The 
grey zone corresponds to values of oj and A for 
which deterministic firing occurs. The values of 
the parameters are / = 1.1, e = 0.05, ipo = 0. 
The initial state is (xo,yo) (courtesy of Pankra- 
tova et al. (0)). 



where e is a fixed small parameter which char- 
acterizes the recovery process, and A, uj are 
respectively amplitude and frequency of the 
external forcing. In the absence of an exter- 
nal driving force there is only one stationary 
state given by 



-I 
-I 



(12) 
(13) 



ex 



;ii) 



The first variable, x, represents the mem- 
brane potential which is characterized by a 
fast dynamics: after x crosses the threshold 
value (x = 0), rapidly it takes on values be- 
low the threshold and a period occurs during 
which no firing is possible. This condition 
is connected with the dynamics of the second 
variable, y, characterized by a slow dynamics. 
Depending on the value of I the stationary 
point is unstable with stable periodic solution 
(|/| < 1) or stable with all the trajectories 
converging on it (attractor) (|/| > 1) (fiol ). 
Here we set / = 1.1 and we take (xo,i/o) as 
initial state for the evolution of the FHN sys- 
tem. In this condition a time evolution oc- 
curs if the system is driven out of equilibrium. 
First, we consider the system in deterministic 
dynamical regime, that is, when firing events 
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FIG. 3 (a) The membrane potential x exhibits 
a spike for A = 0.5 and to = 0.02 (deterministic 
firing region), (b) For A = 0.5, to = 0.01 no spike 
appears (deterministic no- firing region) : the val- 
ues of the parameters don't allow the system to 
fire. Parameter values and initial conditions are 
the same of Fig. [2j 
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FIG. 4 (a) and (b) Stochastic dynamics and 
corresponding trajectory in the space (x, y) for 
A = 0.5, u! = 0.01: because of the coopera- 
tive action of external signal and fluctuations a 
stochastic evolution appears that allows the sys- 
tem to fire out of the deterministic firing region. 
Parameter values and initial conditions are the 
same of Fig. [2j 



occur in the absence of noise (deterministic 
firing region). This regime takes place when 
values of amplitude and frequency are chosen 
inside the gray zone of the parameter plane 
(see Fig. [2]) and corresponds to the situation 
shown in Fig. [3] (a)). Otherwise, when the 
neuron shows no excitability in the absence 
of noise, the system is in the deterministic 
no- firing region (see Fig. [3] (b)). This condi- 
tion occurs when the values of A and lj are 
chosen out of the gray zone of Fig. [2l 



V. THE STOCHASTIC FHN MODEL 

Physical and biological systems are af- 
fected by the presence of continuous random 
perturbations, due to fluctuations of environ- 
mental parameters such as temperature and 
natural resources, which contribute to modify 
the system dynamics (li). These fluctuations 
can be modeled by inserting a noise term in 
the deterministic equations. In this work we 
consider the neuron as an open system whose 
dynamics is affected by the presence both of 
periodic and random variations of environ- 
mental parameters. Therefore we modify the 
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deterministic FHN model by adding a col- 
ored noise term. In particular we consider 
two different cases: i) the membrane poten- 
tial is subject to a noisy dynamics (Case I); ii) 
the refractory variable, which describes the 
recovery properties of the neuron, is exposed 
to fluctuations (Case II). From Eqs. ffTOj) . ffTTj) 
we get 



Case I 

X 

y 



X - — - y + Asin{ujt + ipo) + C 



e(x + J), 



(14) 



Case II 



X 



X = X y + Asin{uj t + (po) 

3 

y = e(x + /) + C; 



(15) 



where (i{t) {i = x,y) are self-correlated 
noises described by Ornstein-Uhlenbeck pro- 
cesses ( j20l ) 

^ = --ut) + -m. (16) 

at T T 



noise enhanced stability (NES). In particu- 
lar the presence of a noise source is responsi- 
ble for modifications in the resonant activa- 
tion phenomenon and causes noise enhanced 
stability to appear in the minimum. NES 
effect and modifications in RA, observed in 
the deterministic firing region (grey zone in 
Fig. 2), are connected with the appearance of 
neuron excitability even if the parameters A 
and u take values out of this region. In fact, 
by setting A = 0.5, u = 0.01 with r ^ 
(white noise) and = 0.005 (Case I), the 
FHN system can fire also for parameter val- 
ues chosen in the deterministic no-firing re- 
gion (see Fig. Hj). Because of the cooperative 
action of periodic force and noise, a stochas- 
tic evolution appears that allows the system 
to fire out of the deterministic firing region. 
However, a self-correlated noise represents a 
model more suitable to represent the random 
fluctuations in real systems, where the noise 
spectrum is characterized by the presence of a 
cut-off. Therefore, in this paper we study the 
system response, in the presence of a periodic 
driving signal, by using the colored Gaussian 
noise given by Eq. (fT6!) (12 ll ). 



Here r is the correlation time of the noise, 
and ^i(t) are statistically independent Gaus- 
sian white noises with zero mean and corre- 
lation function 

< Ut)^j{t') >= v^v^f^ij 6{t-t') = X 

(17) 

The correlation function of the process given 
by Eq. is 



with 



2t 



Le-\^-''\/-6,, (18) 



limr^o{(i{t)(j(t')) = y/o~iy/a]6(t - t')6ij. 

(19) 

In a previous work the system response has 
been studied, under the influence of a peri- 
odic driving signal, function of the white 
noise intensity (fiol ). The authors showed the 
occurrence of resonant activation (RA) and 



VI. RESULTS 

In order to investigate the dynamics of 
I^N system we analyze the behaviour of the 
mean response time (MRT) of the neuron 
when periodical stimulus and colored noise 
are present. Therefore we don't consider the 
periodicity of the signal. MRT is defined as 
< T > = 1/A^ Eili Ti- Here Ti is the first re- 
sponse time (that is, the time for which the 
first spike occurs) of the i*^ realization and 
N is the total number of realizations. In or- 
der to calculate MRT we solve Eqs. f|T^ [T5|) 
by performing numerical simulations, with 

= 5000 for case I (the membrane poten- 
tial is subject to fluctuations) and A^ = 15000 
for case II (the refractory variable is noisy). 
For all realizations the amplitude value of the 
external driving force and the initial condi- 
tions are A = 0.5 and {xo,yo), respectively. 



6 



Case I 




10-3 JQ-2 lQ-1 jqO JQl 10-3 10-2 10-1 lo" 10' 

CO CO 

FIG. 5 (a) and (c): For low noise intensities a slight displacement of the RA minimum is observed 
as the correlation time r increases (weak suppression of the noise effects) . (b) and (d) : For higher 
noise intensities the minimum of RA is strongly affected by the value of r. Both in cases I and II 
the RA minimum almost disappears for white noise, while it is more pronounced for high values 
of r (strong suppression of the noise effects) . Parameter values and initial conditions are the same 

of Fig. m 



We consider a spike occurred when x gets 
over the threshold value Xthr = 0. We ob- 
serve Resonant Activation (RA) and Noise 
Enhanced Stability (NES) phenomena both 
in cases I and II. 



A. Resonant Activation 

In Fig. [5] we report MRT as a function 
of uj for different values both of noise inten- 
sities (Tx-j CTy and correlation time r. The 
value of the RA minimum is affected by the 
noise intensity. For a low level of noise (see 
Fig. [S](a), (c)) this minimum is well pro- 
nounced, without any significative modifi- 
cations occur as the correlation time r in- 
creases. We name this weak suppression of 
the noise effects. However, at higher noise 
intensities, ax = 0.5 and ay = 0.5, (see 



Fig. [5]^b), (d)), the correlation time becomes 
more relevant, so that the value of the RA 
minimum depends strongly on r. In partic- 
ular, as T increases, the minimum is closer 
to the values of the deterministic regime. In 
this case we say that a strong suppression 
of the noise effects occurs. In particular, in 
case I we find a nonmonotonic behaviour of 
MRT function of r in the frequency range 
0.1 < LJ < 1.5. 

B. Noise Enhanced Stability 

By comparing panels (a) and (b) in Fig. \5\ 
for r = (white noise) and r = 5, 10, we ob- 
serve, in case I (membrane potential is sub- 
ject to fluctuations), an enhancement of the 
MRT due to the noise, that is, the depth of 
RA minimum is reduced as the noise inten- 
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sity increases. This enhancement of MRT 
indicates that higher levels of noise cause a 
"response delay" and reduce the neural ef- 
ficiency. Extensively we consider this be- 
haviour as NES effect. A comparison between 
panels (c) and (d) of Fig. O shows that the 
same phenomenon is present also in case II 
(refractory variable is noisy). In view of 
a better understanding of this effect, we an- 
alyze further the behavior of the mean re- 
sponse time. In particular, we calculate the 
MRT as a function of the noise intensity for 
different values of the correlation time. The 
results are shown in Figs. El El We note that 
a weakly correlated noise source affects sig- 
nificatively the efficiency of the neuronal re- 
sponse. In particular, for r = 0.1, 1 we ob- 
serve that, both in cases I and II, the depth 
of RA minimum decreases for high levels of 
noise, which is an enhancement of MRT due 
to the noise (NES effect) (see Figs. and 
[7^,b). Conversely, in the presence of strongly 
correlated noise, no lack of efficiency is ob- 
served around the RA minimum, which main- 
tains its depth almost unchanged as noise in- 
tensity increases (see Figs. [6t,d and[71i). This 
behavior can be evidenced by taking into ac- 
count the values of MRT in the RA minimum. 
From Figs. M and [7] we note that the posi- 
tion of the minimum, for the lowest values of 
the correlation time and the noise intensity 
(r = 0.1 and a = 0.005), is given by ~ 1 for 
case I and u; ~ 0.7 for case II. Therefore we 
calculate MRT function of the noise in- 
tensity for different values of r, setting uj = 1 
in case I and a; = 0.7 in case II. The results 
are reported in Fig. [HI From the inspection 
of the figure we note that, in both cases, the 
strong dependence of the MRT on the noise 
intensity is suppressed as correlation time in- 
creases: the effect of the noise in the RA min- 
imum (region of highest efficiency) appears to 
be smaller as r increases. More precisely we 
note that, in case I, MRT shows in the mini- 
mum a nonmonotonic behavior function 
of a^'- starting from = 0.005 for interme- 
diate levels of noise intensity an increase of 



MRT appears. For higher noise intensities 
a decrease of MRT is observed (see Fig. [8^1). 
On the other side, for case II, a monotonic in- 
crease of MRT is observed clS 8b function of the 
noise intensity cTy. However, in both cases, 
for high values of the correlation time, MRT 
tends to be constant as the noise intensity in- 
creases, in the range of the noise intensity val- 
ues investigated. The modifications induced 
by the noise in the mean response time are 
strongly reduced and they almost disappears 
for r = 10 (suppression of the noise effects). 
Finally we note, in case I, a nonmonotonic 
behaviour of MRT as a function of the cor- 
relation time r. This behaviour appears for 
noise intensity values greater than ax = 0.07, 
corresponding to the maximum of the curve 
for r = 0.1. In Fig. [9l we report MRT vs r 
both for case I and II with a = 0.5. 



C. Colored noise: rescaling effect 

Let's consider two colored Gaussian noise 
sources, 51 and 5*2, with the same intensity 
0", characterized by the correlation times ti 
and T2 respectively, with ti < T2. The be- 
haviour shown in Fig. [HI indicates that the 
presence of a self-correlation causes a reduc- 
tion of the noise "efficacy". This suggests 
that, for frequency values around the RA 
minimum, the noise effect on a FHN system is 
more intensive when one uses 5*1, that is the 
noise source characterized by the smaller cor- 
relation time. Remembering that the white 
noise is obtained by the colored noise in the 
limit r — > (see Eq. (fT^ ). one expects that, 
for frequency values around the RA mini- 
mum, a white noise source with intensity a 
influences the neural dynamics more than any 
correlated noise source with the same inten- 
sity a. 

From Eq. fll9p we note that for |t — 1'| < r 
we can approximate e^'*^*''/^ ~ 1. So that 
from Eq. f[T8l) we get 

mm')) ^ (20) 
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FIG. 6 Case I. For weakly correlated noise (r = 0.1) MRT shows a nonmonotonic behaviour as 
a function of the noise intensity (see panel (a)). For higher values of the correlation time (r = 1) 
this effect is reduced (see panel (b)). In the presence of strongly correlation (r = 5, r = 10) 
this phenomenon disappears and the depth of RA minimum becomes almost independent on the 
noise intensity (suppression of the noise effect)(see panels (c), (d)). Parameter values and initial 
conditions are the same of Fig. [2j 



This result could be interpreted as a rescaling 
effect: the correlation time seems to reduce 
the noise intensity of the colored noise by a 
factor equal to 2r. This suggests that, in the 
RA minimum, the effect of a colored noise 
source with intensity acoiored and correlation 
time r should be, for \t — t'\/T ~ 0, the same 
obtained by using a white noise source with 
intensity given by 

^colored /n-t \ 

(^white = • [^^) 

It 

Therefore we define the "effective" intensity 
of a colored noise according to Eq. fl^Tl) , when 
MRT is of the same order of magnitude of r 
or less. In order to verify this rescaling effect 
we consider again the curve of Fig. [Sji ob- 
tained for Gy = 0.5 and r = 5 and we obtain 
the corresponding "effective" noise intensity 
(^white = {cr colored) I'^r = 0.05. Aftcrwards, by 



using this value we calculate the MRT as a 
function of u. The results have been reported 
in Fig. [TUl In the figure we note that around 
the RA minimum, that is for MRT ~ r and 
^~\t~t'\/T ^ ^Yie influence of the colored 
noise in the FHN dynamics depends on the 
"effective" noise intensity (Twhue = 0.05. This 
means that, in the presence of colored noise 
and for suitable values of the correlation time 
r, the actual noise intensity perceived by the 
FHN system is reduced by the rescaling fac- 
tor l/(2r). We recall that the white noise 
is a theoretical assumption to describe noise 
whose band width is very large. In fact, in 
real systems the fluctuations are connected 
with colored noise sources. Therefore, in or- 
der to evaluate the response of a real neu- 
ron, one has to consider self-correlated noise 
sources: the rescaling phenomenon, as we 
show in this paper, reduces the effect of the 
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FIG. 7 Case II. For weakly correlated noise (r = 0.1) a monotonic behaviour of MRT as a 
function of the noise intensity appears (see panel (a)). As the correlation time increases (r = 1) 
this phenomenon is attenuated (see panel (b)). Finally, for strongly correlation (r = 5, r = 10) 
the depth of RA minimum is almost independent on the noise intensity (suppression of the noise 
effect). Parameter values and initial conditions are the same of Fig. [2j 



noise on the FHN dynamics, modifying in a 
significative way the neuronal dynamics. 



VII. CONCLUSIONS 

After a brief discussion on the Hodgkin- 
Huxley (HH) model, we introduce the 
FitzHugh-Nagumo (FHN) model. Although 
it is a simplified version of the former, the 
latter permits to separate the mathemati- 
cal characteristics of excitation and propaga- 
tion from the electrochemical peculiarities of 
sodium and potassium ion flows. Even if the 
HH model better describes the real dynamics 
of the neuronal response, it allows only to ob- 
serve two-dimensional projections of its four- 
dimensional phase trajectories. Because of 
this, the FHN model is often preferred, pro- 
viding the whole solution directly in a two- 
dimensional phase space. This characteris- 



tic permits to envisage a geometrical inter- 
pretation of important biological phenomena 
that depend on the spike-generating mecha- 
nism which causes the neuronal response to 
external stimuli. Therefore, in this paper we 
analyzed the response of a neuron in the pres- 
ence both of a driving periodical force with 
frequency uj and an additive Gaussian col- 
ored noise, by the FHN neuronal model. In 
our analysis we considered two cases: I) the 
noise source affects the membrane potential, 
II) the noise source influences the recovery 
variable. In these conditions we analyzed 
the mean response time (MRT) as a function 
both of the noise intensity and the correlation 
time. We found that RA and NES phenom- 
ena undergo meaningful modifications due to 
the presence of different values of the correla- 
tion time r. For strongly correlated noise we 
observed suppression of NES and persistence 
of RA (efficiency enhancement of neuronal re- 
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FIG. 8 (a) Case I. Nonmonotonic behaviour 
(with a maximum) of MRT as a function of the 
noise intensity ax, for ui = 1.0: in the presence of 
colored noise this effect disappears as correlation 
time increases, (b) Case II. Monotonic increase 
of MRT as a function of the noise intensity ay, 
for Lo = 0.7: this behaviour is suppressed for 
strongly correlated noise. Parameter values and 
initial conditions are the same of Fig. [2l 



sponse) as correlation time r increases. This 
indicates that the RA minimum correspond- 
ing to a certain value of u is preserved when 
the noise source is characterized by large val- 
ues of the correlation time. Conversely the 
enhancement of the MRT, which indicates a 
significant role of the noise, tends to van- 
ish. The reduction of the noise effects in 
the presence of strongly correlated noise in- 
dicates a rescaling effect of the noise self- 
correlation. We investigated this aspect for 
a colored noise source whose intensity and 
correlation time are given by a and r respec- 




FIG. 9 (a) Case I. Nonmonotonic behaviour of 
the MRT as a function of the correlation time r, 
for uj = 1.0 and ax = 0.5. (b) Case II. Decrease 
of the MRT as a function of the correlation time 



T, for CO = 0.7 and a. 



y 



0.5. Parameter values 



and initial conditions are the same of Fig. [2j 



tively. We found that, for MRT ~ r, an 
"effective" noise intensity exists. This im- 
plies that the effect of the self-correlated noise 
can be reproduced, through a rescaling pro- 
cedure, by using a white noise source whose 
amplitude is given by ay^hue = {cr colored) /'2^r. 
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FIG. 10 Case II. The correlation time causes a 
rescaling effect according to au^hUe = cr colored f^r. 
Around the minimum the curve obtained for 
(^colored = 0.5 and correlation time r = 5 
(black dash line) overlaps that obtained by us- 
ing a white noise source with intensity awhUe = 
colored /"^T = 0.05 (gray solid line). Parame- 
ter values and initial conditions are the same of 
Fig. [21 



mazione (PON 2000-2006). More informa- 
tion is available at http:/ /ww w.pi2s2.it, and 
[http: //www.consorzio-cometa.it 



References 

[1] See the special section on "Complex Sys- 
tems", Science 284, (1999) 79-107; O. 
N. Bjornstad and B. T. Grenfell, SCI- 
ENCE 293, (2001) 638-643; S. Ciuchi, F. de 
Pasquale and B. Spagnolo, Phys. Rev. E 53, 
(1996) 706-716; M. Scheffer, S. R. Carpen- 
ter, J. A. Foley, C. Folke, B. Walker, Na- 
ture 413, (2001) 591-596; B. Spagnolo, D. 
Valenti, A. Fiasconaro, Math. Biosc. and 
Engineering 1, (2004) 185-211. 

[2] C. R. Doering and J. C. Gadoua, Phys. Rev. 
Lett. 69, (1992) 2318-2321; R. N. Man- 
tegna and B. Spagnolo, Phys. Rev. Lett. 84, 
(2000) 3025-3028; P. Majee, G. Goswami, 
B. Chandra Bag, Chem. Phys. Lett. 416, 
(2005) 256-260; R. Gommers, P. Douglas, 
S. Bergamini, M. Goonasekera, P. H. Jones, 



and F. Renzoni, Phys. Rev. Lett. 94, (2005) 
143001(1-4). 

[3] R. N. Mantegna, B. Spagnolo, Phys. Rev. 
Lett. 76, (1996) 563-566; R. Wackerbauer, 
Phys. Rev. E 59, 2872 (1999); A. Mielke, 
Phys. Rev. Lett. 84, (2000) 818-821; B. 
Spagnolo, A. A. Dubkov, and N. V. Agu- 
dov. Acta Phys. Pol. 35, (2004) 1419-1436; 
G. Bonanno, D. Valenti and B. Spagnolo, 
Phys. Rev E 75, (2007) 016106(1-8). 

[4] N. V. Agudov, B. Spagnolo, Phys. Rev. E 
64, (2001) 035102(1-4); A. A. Dubkov, N. 
V. Agudov and B. Spagnolo, Phys. Rev. E 

69, (2004) 061103(1-7). 

[5] E. Lanzara, R. N. Mantegna, B. Spag- 
nolo and R. Zangara, Am. J. Phys. 65, 

(1997) 341-349; L. Gammaitoni, P. Hanggi, 
P. Jung, F. Marchesoni, Rev. Mod. Phys. 

70, (1998) 223-288; Y. Kashimori, H. Fu- 
nakubo, and T. Kambara, Biophys. J. 75, 

(1998) 1700-1711; D. Valenti, A. Fias- 
conaro, B. Spagnolo, Phys. A 331, (2004) 
477486; B. Kosko and S. Mitaim, Phys. 
Rev. E 70, (2004) 031911 (1-10); A. Caruso, 
M. E. Gargano, D. Valenti, A. Fiasconaro, 
B. Spagnolo, Flue. Noise Lett. 5, (2005) 
L349-L355. 

[6] K. Wiesenfeld, F. Moss Nature 373, (1995) 
33-36; F. Chapeau-Blondeau, X. Godivier, 
and N. Chambet Phys. Rev. E 53, (1996) 
1273-1275 ; Special Issue: Advances in neu- 
ral networks research IJCNN'03, B. Kosko, 
S. Mitaim, Neural Networks 16, (2003) 755- 
761. 

[7] A. L. Hodgkin and A. F. Huxley, J. Phys- 
iol. 117, (1952) 500-544; E.V. Pankratova, 
A.V. Polovinkin, E. Mosekilde, Eur. Phys. 
J. B 45, (2005) 391-397; V.N. Belykh and 
E.V. Pankratova, Dynamics and synchro- 
nization of noise perturbed ensembles of pe- 
riodically activated neuron cells, Int. J. Bi- 
furcation and Chaos 18, (2008) in press. 

[8] K. F. Bonhoeffer, J. Gen. Physiol. 32, 
(1948) 69-91; K. F. Bonhoeffer, Naturwis- 
senschaften 40, (1953) 301-311. 

[9] B. van der Pol and J. van der Mark, Arch. 



12 



Neerl. Physiol. 14, (1929) 418-443 . 

[10] R. FitzHugh, Bull. Math. Biophysics 17, 
(1955) 257-278; R. FitzHugh, J. Gen. Phys- 
iol. 43, (1960) 867-896; R. FitzHugh, Bio- 
phys. J. 1, (1961) 445-466. 

[11] J. H. Hale, H. Kocak, Dynamics and Bifur- 
cations (Springer- Verlag, New York 1991). 

[12] C. Rocsoreanu, A. Georgescu, N. Giur- 
giteanu. The FitzHugh- Nagumo Model: Bi- 
furcation and Dynamics (Kluwer Academic 
Pubhshers, Boston 2000). 

[13] J. Nagumo, S. Arimoto, and S. Yoshizawa, 
Proc. IRE 50, (1964) 2061-2070. 

[14] J. Keener, J. Sneyd, Mathematical Physiol- 
ogy (Springer- Verlag, New York 1998). 

[15] X. Pei, K. Bachmann and F. Moss, Phys. 
Lett. A 206, (1995) 61-65. 

[16] A. S. Pikovsky, and J. Kurths, Phys. Rev. 



Lett. 78, (1997) 775-778. 
[17] D. Nozaki, Y. Yamamoto, Phys. Lett. A 

243, (1998) 281-287. 
[18] A. Longtin, D. R. Chialvo, Phys. Rev. Lett. 

81, (1998) 4012-4015. 
[19] E. V. Pankratova, A. V. Polovinkin, B. 

Spagnolo, Phys. Lett. A 344, (2005) 43-50. 
[20] C. W. Gardiner, Handbook of Stochastic 
Methods for Physics, Chemistry and the 
Natural Sciences, (Springer, Berlin 1993). 
[21] D. Valenti, A. Fiasconaro and B. Spagnolo, 
Mod. Prob. Stat. Phys. 2, (2003) 91-100; 
D. Valenti, A. Fiasconaro and B. Spagnolo, 
Flue. Noise Lett. 5, (2005) L337-L342; D. 
Valenti, L. Schimansky-Geier, X. Sailer, B. 
Spagnolo, M. lacomi, A. Phys. Pol. B 38, 
(2007) 1961-1972. 



13 



